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Abstract 

Many recursive functions can be defined elegantly as the unique homomorphisms, between 
two algebras, two coalgebras, or one each, that are induced by some universal property 
of a distinguished structure. Besides the well-known applications in recursive functional 
programming, several basic modes of reasoning about scientific models have been demonstrated 
to admit such an exact meta-theory. Here we explore the potential of coalgebra-algebra 
homomorphism that are not a priori unique, for capturing more loosely specifying patterns 
of scientific modelling. We investigate a pair of dual techniques that leverage (co)monadic 
structure to obtain reasonable genericity even when no universal properties are given. We 
show the general applicability of the approach by discussing a suprisingly broad collection of 
instances from real-world modelling practice. 


1 Introduction 

This paper explores a (co) algebraic framework for homomorphic and recursive reasoning in and 
about scientific models; that is, mathematical structures we think in, loaded with interpretations of 
phenomena from the world we live in. To this end, the paper is structured as follows: the remainder 
of section 1 reviews the relevant concepts and notations of categorial (co)algebra, with an interlude 
in section 1.3 that states the motivation and goal more precisely, once terms have been established. 
Section 2 introduces our dual pair of recursion schemes of interest. Sections 3 and 4 discuss two 
application domains per scheme. These should be regarded as a set of four related short papers, 
and are mostly self-contained with respect to discussion and related work. A general conclusion is 
difficult because of the wide variety of scope, and has to be omitted due to space constraints. We 
trust the application subsections to speak for themselves. Proofs and illustrations are relegated to 
the appendix. 


1.1 (Co)Algebras of a Functor 

Let F be an (endo)functor on a category C , in all of the following, but without loss of generality, 
the category of sets. F-algebras are structures (X,f : FX —> X ), where X is called the carrier 
and / the operation. U-coalgebras are the dual structures {X,f : X —> FX). Homomorphisms 
are morphisms that make a square commute. They can be defined between two algebras or two 
coalgebras, or between one of each: 
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The former two cases of pure (co)algebra homomorphisms have been studied quite extensively in 
the frameworks of universal (co)algebra. The latter two mixed cases have only relatively recently 
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received significant attention. The last case, of homomorphisms from an algebra to a coalgebra, 
remains rather obscure. By contrast, homomorphisms from a coalgebra to an algebra, called 
coalgebra-algebra homomorphisms and in the following abbreviated as ca-homomorphisms, are 
a very general and expressive framework for divide&conquer schemes of recursion. The common 
intuition (paraphrased from [5]) is the following: ‘In order to solve ( h ) a computation problem 
on complex inputs (X), decompose (/) them into a collection of subproblems (XA), solve these 
independently ( Fh ) to obtain a collection of subresults (XY), and compose ( g ) those to form the 
final result (Y).’ 

The composition of two F-(co) algebra homomorphisms is again an X-(co)algebra homomorphism, 
respectively. The composition of an X-algebra homomorphism with an X-ca-homomorphism, or du¬ 
ally an X-ca-homomorphism with an X-coalgebra homomorphism, is again an X-ca-homomorphism. 
The following proposition about “new ca-homomorphisms from old ones” is less obvious, but useful 
for the subsequent discussion. 

Lemma 1.1. Let functor X preserve (co)products. Then the (co)tuplings of ca-homomorphisms 
are the ca-homomorphisms associated with (co)product-structured (co)algebras, respectively , up to 
unique isomorphism. 

1.2 Distinguished (Co)Algebras 

The X-algebras and X-algebra homomorphisms form a category Alg(X). Of particular interest 
are initial X-algebras (/^X, in). 1 For any X-algebra (A, /) there is a unique homomorphism 
fl/D : (//X, in) — > (X, /), called a catamorphism. By Lambek’s lemma, the operation (in) is a 
bijection. Dually, the X-coalgebras and X-coalgebra homomorphisms form a category Coalg(X). 
Of particular interest are final X-coalgebras (t'X, out). For any X-coalgebra (A', /) there is a unique 
homomorphism [/] : (A, /) —► (i'X, out), called an anamorphism. The operation (out) is a bijection. 
For initial algebras and final coalgebras of a variety of well-behaved functors, see [16]. 

For ca-homomorphisms, the situation is more complex, because they do not compose among 
themselves. A X-algebra (Y, g) is called corecursive [6] if and only if it plays a role analogous to a final 
X-coalgebra: For any X-coalgebra (A, /) there is a unique X-ca-homomorphism h : (A, /) —► (Y, g). 
Dually, a X-coalgebra (A, /) is called recursive [5, 12] if and only if it plays a role analogous to an 
initial X-algebra: For any X-algebra (Y, g) there is a unique X-ca-homomorphism h : (A, /) —> (Y, g). 
Trivial examples are obtained by Lambek’s lemma: (/xX, in -1 ) is a recursive X-coalgebra; and 
dually (uF, out -1 ) is a corecursive X-algebra. 

1.3 Interlude: Vision 

The present paper deals with intermediate situations: Neither are our ca-homomorphisms of interest 
completely arbitrary, with algebra and coalgebra chosen ad hoc at the same level of particularity; 
nor is there an obvious candidate for a general (co)recursive (co)algebra, leaving only the partner to 
be chosen in particular, and the homomorphism induced uniquely. Instead, we investigate situations 
where either the algebra or the coalgebra part is more central, and the possibility and general shape 
of homomorphisms is to be studied, independently from its less fixed partner. 

This line of investigation ties into our overarching research programme on the recursive nature 
of scientific modelling posited by Rosen [14]. We shall identify certain modes of abstract formal 
reasoning over classes of models as (co)algebras, and their rigorous interpretation as “reusable” 
ca-homomorphisms between a fixed general formal language and a variety of particular models. 

In previous work, we have investigated dual, purely algebraic and coalgebraic reasoning modes, 
and identified them as ‘queries of causality’ and ‘representations of behavior’, respectively [8]. 
We have also explored various types of mixed modes, and drawn connections between structural 
operational semantics and cellular automata [18], and between course-of-value iteration and history- 
dependent dynamics [19], respectively. The present work deals with yet more types of mixed mode, 
but ones that come as a dual pair, as in [8]. We consider the duality to be of great philosophical 
interest, although the present modest paper can only give a few first directions. 

1 Note the boldface typography that distinguishes fixpoint operator //, and monad multplication //. 
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A great deal of the desired reusability shall be achieved by imposing a (co) monad structure on 
the functor T. In the following section, we review the relevant concepts and notations very briefly, 
aware that our presentation can only serve as a glossary, not an introduction. 

1.4 (Co)Monads 

A monad T = (T , 77 , p) is an endofunctor T together with two natural transformations 77 : 1 => T 
and p : TT =>■ T, called unit and multiplication , respectively, such that p o r)T = idT = p oTp and 
p o pT = p oTp. A monad gives rise to a Kleisli extension operator that sends any morphism 
/ : A —> TY to a unique /* : TX —> TY, such that 77 * = idT, /* o p x = f and g* o /* = (< 7 * o /)*; 
namely /* = p Y 0 Tf and / = /* o r) X . 

The Kleisli category K1(T) of a monad has the same objects as the underlying category, but 
homomorphisms / : X K1 ^ T \ Y whenever / : X —> TY. Its identity is id K1 * T ^ = 77 ; composition is 

given by g o K h T ) f = g* o f. Coalgebras of a monadic functor are of Kleisli type / : X —X 
and can be iterated. We write /* ,n for the morphism such that f* ,n = (/*)". 

For many nice endofunctors T, a free monad F* can be constructed as follows [3] , 

F*X = p(X + F ) T*(/ : X -> Y) = <\in Y o (/ + id FF .y)|) 

F* is a functor, and turns the family of initial algebras into a natural transformation in : Id+TT* => 
F*. Besides the usual natural transformations that go with a monad, define two additional ones, 
r : FF* => F* and k : F => F* (adapted from [4]), as well as Kleisli extension, simultaneously as: 

77 = in o i 1 t = in o t 2 p = (idT*)* k = t o T ?7 (/ : X —> F*Y)* = (] [/, tv] ) 

Lemma 1.2. The extra transformations obey the law r = p o kF* . 

Dually, a comonad is an endofunctor D together with two natural transformations e : D => 1 and 
v : D =>• DD, called counit and comultiplication, respectively, such that eD o v = idT = Ds o v and 
uD o v = Dv o v. A comonad gives rise to a co-Kleisli extension operator that sends any morphism 
/ : DX —> Y to a unique f° : DX —> DY , such that e° = idT, e Y o f° = / and g° 0 f° = (g 0 /°)°; 
namely f° = Df o v x and f = e Y o f°. 

The co-Kleisli category Cl(O) of a comonad has the same objects as the underlying category, 

but homomorphisms / : X y whenever / : TA' —> Y\ Its identity is id cl *- D ^ = e; composition 
is given by g o cl ( D ) / = g o f°. 

A cofree comonad T°° can often be constructed from an endofunctor F as follows. 

F°°X = u{X x T) F°°(f : F°°X ->Y) = [((/ x id FF ~*) o out*)] 

T°° is a functor, and turns the family of final coalgebras into a natural transformation out : T°° => 
Id + TT°°. Besides the usual natural transformations that go with a comonad, define two additional 
ones, 6 : T°° => TT°° and x : => T (dually extrapolated from [4]), as well as co-Kleisli 

extension, simultaneously as: 

£ = 7 Tioout 0 = 7 t 2 o out v = (idT 00 ) 0 x = Te o 0 (/ : T°°A -4 F)° = [(/, 0 X )J 

Lemma 1.3. The extra transformations obey the law 0 = xF°° 0 v - 

1.5 (Co)Algebras of a (Co)Monad 

An algebra (A',/) of a functor T is also an algebra of the monad T = (T,rj,p), also called 
an Eilenberg-Moore algebra, if and only if / is compatible with the monad operations, that is 
/ 0 Vx = id.Y and /o p x = f oTf. The algebras of a monad Alg(T) form a full subcategory of the 
algebras Alg(T) of the underlying functor. The monad multiplication is a distinguished algebra 
operation: 

Lemma 1.4. The T -algebras of the form (TX,p x ) are “locally” weakly initial: for every T -algebra 
(A,/), there is a canonical T-algebra homomorphism f : ( TX,p x ) —» (A,/). Furthermore , any 
morphism h : X Y gives rise to a T-algebra homomorphism Th : ( TX : p x ) —> (TY,p Y ). 
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The algebras of a functor F and the algebras of the free monad F* = (F* ,rj, /i) are in one-to-one 
correspondence: Each F*-algebra (A,/) is of the form / = ef = (] [idx, e] [) , generated uniquely by 
the F-algebra ( X , e) where e = / b = / o kx ■ Thus we can reformulate Lemma 1.2. 

Lemma 1.5. The free monad extra transformations obey the law r = \j?. 

Dually, a coalgebra (A', /) of a functor D is also a coalgebra of the comonad D = ( D , e, i/), if 
and only if £x o / = idx and vx ° f = Df 0 f ■ The coalgebras of a comonad Coalg(O) form a full 
subcategory of the coalgebras Coalg(D) of the underlying functor. 

Lemma 1.6. The D -coalgebras of the form (M, g\-) are “locally ” weakly final: for every O- 
coalgebra (A,/), there is the canonical D-coalgebra homomorphism f : (A,/) —> {DX,vx)- Fur¬ 
thermore, any morphism h : X —» Y gives rise to a D-coalgebra homomorphism Dh : (DA, vx) —> 
{DY, v Y ). 

The coalgebras of a functor F and the coalgebras of the cofree comonad F°° = (F°°,e,v) are in 
one-to-one correspondence: Each F°°-coalgebra (A,/) is of the form / = ef = [ (id,Y, e) ], generated 
uniquely by the F-coalgebra (A, e) where e = f b = \x 0 /• 

Lemma 1.7. The cofree comonad extra transformations obey the law 6 = . 


2 (Co)Kleisli (Co)lnduction 

We shall demonstrate that well-behaved classes of interesting ca-homomorphisms arise from a pair of 
dual recursion schemes that leverage (co)monadic structure. The monadic version has been studied 
recently in the context of trace semantics [4], Although its dual comonadic twin is straightforward, 
we are not aware of previous uses; in fact we shall argue why its characteristic properties, albeit 
essentially useful for our purpose here, are undesirable in many circumstances. 

2.1 Kleisli Coinduction 

Definition 2.1 (Kleisli Coinduction). Let T = (T, 77 ,/x) be a monad. Let (A, e) be a T-coalgebra. 

A morphism e' : X —» TY is called a Kleisli-coinductive solution of equation e, if and only if 
e t = e t 0 K1(T) e _ 

Lemma 2.2 (Characterization). Kleisli-coinductive solutions e' : X — > TY are precisely the 
T-ca-homomorphisms to the multiplicative T -algebra (TY^fiy)- 

A Kleisli-coinductive solution e' : X —> TY of e : A — > TX appears as a morphism of type 
A —► Y in the Kleisli category. As usual for corecursive functions, its codomain Y is not determined 
by e at all. 

Lemma 2.3 (Substitution). Let (A, e) be a T-coalgebra and h : Y Z be any morphism. If 
e' : A' —> TY is a Kleisli-coinductive solution, then so is Th o e' : X —» TZ. 

Our motivation for studying this particular corecursion scheme is that each Kleisli-coinductive 
solution extends to many ca-homomorphisms in a regular way, such that its common properties 
can be studied regardless of the target algebra. 

Theorem 2.4 (Universality). Given any T-coalgebra (A, e) and Kleisli-coinductive solution e' : 
A —> TY, there is a canonical family of extensions to T-ca-homomorphisms (/ o C) into all 
T -algebras ( Y, f ). 

Corollary 2.5 (Recursivity). If a T-coalgebra (A, e) is recursive, then it has a unique Kleisli- 
coinductive solution, whose extension is the unique ca-homomorphism. 

Conversely, if a T-coalgebra has a unique Kleisli-coinductive solution, then it is “morally 
recursive”: there is a unique ca-homomorphism into any T-algebra that is canonical in the sense of 
Lemma 2.2. Note that T-algebras which are incompatible with the monadic structure are generally 
not covered. 
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On the other hand, Kleisli-coinductive solution can be non-unique, or fail to exist at all. 
Hence one should regard the equation coalgebra as a specification rather than definition of ca- 
homomorphisms. Such a specification can be contradictory or loose; the following proposition shows 
that it may also be tautological. 

Lemma 2.6 (Unit). Let T = (T, fy,/i) be a monad. For a T-coalgebra of the form (X,rjx), all 
morphisms of type X —> TY are Kleisli-coinductive solutions. 

While a globally tautological specification is evidently of little practical value, we shall make 
good use of pointwise local tautology in the application section 3.1. 

For the special case of a free monad F*, we have an asymmetric situation: The E*-algebras have 
a distinguished subclass of interest, namely the F*-algebras which are equivalent to the simpler 
class of E-algebras. But F* -coalgebras generally obey no such constraint of expressivity. If one is 
imposed deliberately, then a reduced characterization of Kleisli-coinductive solutions in terms of 
E-ca-homomorphisms can be given. 

Definition 2.7 (Basic Equation). Let F* = (F* ,r],p) be a free monad over endofunctor F. An 
E*-coalgebra (A, e) is called basic , if and only if it factors as e = Kx o eo- Note that (X, eo) is an 
E-coalgebra. 

Lemma 2.8 (Demonadization). F*-ca-homomorphisms from basic F*-coalgebras (X, Kx ° eo) to 
F *-algebras (Y, /) are precisely the F-ca-homomorphisms from (X, eo) to (Y,f b ). 

Corollary 2.9. The Kleisli-coinductive solutions e' : X —> F*Y of basic F*-coalgebras (X, kx °eo) 
are precisely the F-ca-homomorphisms from (A', eo) into the F-algebra (E*Y, 7y). 

2.2 Co-Kleisli Induction 

All of the preceding statements can be dualized. 

Definition 2.10 (Co-Kleisli Induction). Let D = (D, e, v) be a comonad. Let (Y, k) be a E-algebra. 
A morphism : DX —» Y is called a co-Kleisli-inductive solution of coequation k, if and only if 
fct = ko c m k t. 

Lemma 2.11 (Characterization). Co-Kleisli-inductive solutions kf : DX —> Y are precisely the 
ca-homomorphisms from the D-coalgebra (DXyx). 

Lemma 2.12 (Substitution). Let ( Y,k ) be a D-algebra and h : W —>• X be any morphism. If 
kf : DX —> Y is a co-Kleisli-inductive solution, then so is kf o Dh : DW —> Y. 

Theorem 2.13 (Universality). Given any D-algebra (Y, k) and co-Kleisli-inductive solution k$ : 
DX —> Y, there are canonical extensions to ca-homomorphisms (k$ o /) from all D -coalgebras 
(*,/)• 

Corollary 2.14 (Corecursivity). If a D-algebra (Y, k) is corecursive, then it has a unique co- 
Kleisli-inductive solution, whose extension is the unique ca-homomorphism. 

Lemma 2.15 (Counit). Let D = (D,e,v) be a comonad. For a D-algebra of the form (Y, £y), all 
morphisms of type DX —> Y are co-Kleisli-inductive solutions. 

Definition 2.16 (Cobasic Coequation). Let F°° = (F°°,e, v) be a cofree comonad over endofunctor 
F. An F°°-algebra (Y k) is called cobasic, if and only if it factors as k = fco ° Xy- Note that 
(Y, kg) is an E-algebra. 

Lemma 2.17 (Decomonadization). F°°-ca-homomorphisms from F°°-coalgebras (X,/) to cobasic 
F°°-aZ</e&ras (Y, ko o \y) are precisely the F-ca-homomorphisms from (X, f b ) to (Y, ko). 

Corollary 2.18. Co-Kleisli-inductive solutions kf : F°°X —> Y of cobasic F°°-algebras are precisely 
the F-ca-homomorphisms from the F-coalgebra (F°°X,dx)- 
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Our examples in section 4 build on cofree comonads, which can be understood to create spaces 
of non-well-founded node-labelled trees. In this context, induction seems like a thing out of hell: 
co-Kleisli-inductive solutions are shamelessly allowed to depend on infinite regresses, as well as in 
vicious circles (via e) on their own results. But, as we have pointed out, we are happy to sacrifice 
universal existence and uniqueness of solutions; the scheme is still quite useful as a formal framework 
for organizating existing modes of reasoning, whose potential unsoundness is either of no practical 
concern, or resolved in domain-specific ways. 


3 Applications of Kleisli Coinduction 

3.1 Dynamical Systems 

Choose a monoid, written additively as (A, 0, +). It may or may not be Abelian, or a group. The 
(left) product with a monoid, T = Ax Id, is a monad with: 

r](x) = (0, x) fi(t, (u , a;)) = (t + u, x) f*(t, x) = (t + u, y ) where f(x) = (u , y) 

The algebras of this monad are exactly the (left) monoid actions; (A, / : A x I -> A) where: 

/(0, x) = x f(t + u,x) = f(t,f(u,x)) 

Monoid actions are the most general formal framework for dynamical systems. The monoid 
is understood as the structure of durations of time. Default candidates are non-negative/all 
integers/reals, which model irreversible/reversible discrete/continuous time, respectively. The 
language that arises from T deals exclusively with the passing of time: In covariant positions, such 
as the codomain of morphisms in the Kleisli category, (t,x) £ TX reads as “after t time units, 
x"\ in contravariant positions, such as the domain of T-algebras, (t,x) £ TX reads as “f time 
units after x”. The Kleisli category is the category of sets and timed functions, where f(x) = (t , y) 
consequently reads as “/ maps x after t time units to y”. Kleisli composition is sequential in time; 
delays are accumulated. 

Example 3.1 (Harmonic Oscillator). Consider the textbook example of a harmonic oscillator, a 
point mass m moving frictionlessly along a straight line, acted on by a restoring force proportional, 
with coefficient fc, to its displacement x. In a Newton-style modelling approach, this system is 
specified by a second-order linear differential equation x + ^x = 0, which simply states that all 
acceleration is due to the restoring force. 

The appropriate analytic model is a dynamical system with A = ffi. and state vectors ( x , x) £ M 2 . 


f(t,{x,x)) 


f cos ut 
wsin cot 


CO 


-1 


sin cot 
cos cot 



where co = 



The ca-homomorphism diagram for this monad takes the equation coalgebra (A, e) as a set A of 
symbolic states, together with a map e that gives each state x £ X a backwards-looking specification 
(t,x'), understood as “t time units after state x' v . The ca-homomorphisms h : (A, e) —X (Y, f) are 
then the instantiations consistent with the forwards-looking dynamics / on a concrete state space 
Y. Note that partial specifications are included automatically: A state specified tautologically as 
zero time units after itself is consistent with any instantiation, by pointwise Lemma 2.6. 

By unfolding the definition, the Kleisli-coinductive solutions of a T-coalgebra (A, e) are exactly 
the morphisms of type e* : A — > TY such that e(x) = ( t, x') and e'(x’) = ( u , y) implies e^(t + u,y). 
They have a surprisingly rich structure, and admit necessary and/or sufficient conditions, under a 
variety of mild assumptions concerning the time-likeness of A. As outlined before, our framework 
allows us to study the existence and parameters of solutions, independently of interpretation in a 
concrete dynamical system. 

Definition 3.2. Associate with a T-coalgebra (A, e) several relations: 

• (~ e ) C A x A is the reflexive-symmetric-transitive closure of TT 2 o e, an equivalence. 

• (Ae) C A x A is the reflexive-transitive closure of tt 2 o e, a preorder. 
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— Partition into complementary subrelations (V e ) = (V e ) \ (^ e ), a strict partial order, and 
(— e ) = (^ e ) H (^ e ) C (~ e ), an equivalence, as usual. 

• (~^ e ) CXx TA is the Kleisli analog of (V e ), namely (-w e ) = U^Lo e *’ n - 
We call (~-> e ) consistent, if and only if x (t, x) implies t = 0. 

These relations allow us to state some basic properties of Kleisli-coinductive solutions. 

Lemma 3.3. Let (A, e) be a T-coalgebra. Then it splits into a coproduct of the (~ e )- equivalence 
classes. Furthermore, every Kleisli-coinductive solution e' : X —> TY is locally constant in its second 
component on each (~ e )-equivalence class: If x ~ e x', then e'(x ) = (u,y) and e'(x') = ( u',y') 
imply y = y'. 

Lemma 3.4. If (-w e ) is consistent, then for each pair x > e x' there is a unique t such that 
x ( t,x'). 

Lemma 3.5. Let e' : X — » TY be a Kleisli-coinductive solution of (A, e). If x ( t,xand 
e\x') = ( u,y ), then e\x) = ( t + u,y ). In particular, x (0, cc 7 ) implies e\x) = e'(x'). 

We say that the monoid (A,0, +) has right cancellation, if and only if t + u = u implies t = 0 
for all t,u £ A. For cancellative monoids, a reasonable assumption for standard models of time, we 
can give both a necessary and a sufficient condition on the solvability of a T-coalgebra. 

Theorem 3.6. Assume the monoid (A,0,+) has right cancellation. Then a T-coalgebra (A, e) has 
Kleisli-coinductive solutions, only if (~-> e ) Is consistent. 

Proof. Assume, for contradiction, that x ~> e (t,x) where t / 0. Let (u,y) = e^(x). Then also 
e^(x) = (t + u,y), by Lemma 3.5; hence t + u = u, and finally by cancellation t = 0. □ 

Theorem 3.7. Assume the monoid (A,0,+) has right cancellation. Then a T-coalgebra (A', e) has 
Kleisli-coinductive solutions, if (~-> e ) is consistent, and (>- e ) is well-founded. 

Proof (Sketch). By well-founded induction and Lemma 3.5. □ 

For reversible time, the conditions can be made tight. 

Theorem 3.8. Assume (A,0,+) is a group. Then a T-coalgebra (A', e) has Kleisli-coinductive 
solutions, if (and only if) (-w e ) is consistent. 

Proof (Sketch). Componentwise by Lemma 3.3, then either by reduction to Theorem 3.7 or else by 
inversion of Lemma 3.4. □ 

So far, we have described the solution spaces of T-coalgebras formally, but given no explication 
of their meaning in the modelling context. We shall now demonstrate that the T-ca-homomorphisms 
are a very abstract and general account of time series: discrete samples of state snapshots over the 
dynamic evolution of a system. 

Fix A as either the nonnegative or all reals. Fix a bilaterally infinite sequence of time differences, 
(Si £ R) for all i £ Z. We write (tf) for the corresponding partial sums: 

U = + J2 Si if i > 0 ti = - J2 Si if i < 0 

je[o,i) je[i, o) 

Now consider the following T-coalgebras on integer intervals: 


A-| = (—oo, + 00 ) 


ei(i + 1) = (Si,i) 

(i £ Z) 

A 2 = (— 00 ,0] 


e 2 (i + 1) = (Si,i) 

(*< o) 

A 3 = [ 0 , + 00 ) 

e 3 ( 0 ) = ( 0 , 0 ) 

e 3 (i + 1 ) = (8i,i) 

( 0 <<) 

A 4 = [0,n] 

e 4 ( 0 ) = ( 0 , 0 ) 

e 4 (i + 1) = (Si,i) 

(0 < i < n) 


We shall verify that the induced ca-homomorphisms from e^ into dynamical systems are the different 
shapes of time series: for k = 1,2, 3,4, the bilaterally infinite, the left-infinite, the right-infinite, 
and the finite of length n + 1 , respectively. 
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In each case, there is only one equivalence class of (~ eit ). Hence any Ivleisli-coinductive solution 
e\ : Xk —> TY fixes a single target element yo G Y, and is of the form ej.(i) = (ti + uo,yo). 
Consequently, the ca-homomorphisms h = f o e\ : ( Xk,ek) —t (Y, /) are of the form: h(i) = 
f(ti + u o, yo). That is, they are sequences of elements on the trajectory of reference state yo, spaced 
according to the sampling time sequence (ti) relative to reference time uo ■ Without loss of generality, 
Uq can be made to vanish, by rewriting to h(i) = f(ti,y' 0 ) where y' 0 = f(uo,yo)- Furthermore, if 
(Si) is constant, then we obtain the equidistant time series, with the simpler form h(i) = f(i ■ 5,yo). 

Example 3.9 (Period). For the harmonic oscillator from Example 3.1, the well-known periodic 
behaviour can be captured very concisely and naturally, by stating that the bilaterally infinite, 
equidistant time series specification coalgebra (Ai,ei) with 8 = 2ttuj~ 1 admits only constant 
ca-homomorphisms. 

Example 3.10 (Zeno). As a sidenote, the pathological class of non-equidistant, right-infinite but 
bounded time series have played an important, puzzling role in ancient Greek “scientific modelling”. 
For instance, Zeno’s stadium run can be specified by T-coalgebra (A 3 , 63 ) with Si = 2~^ +1 \ where 
half of the remaining time is consumed at each step. 

This account of discrete time series over continuous dynamical systems may seem a little 
contrived, simply because standard textbook presentations look fairly different, and do not suggest 
the use of a formal framework beyond simple set theory and “index magic”. The following example 
shall serve, among other things, to counter that impression; there the study of ca-homomorphisms 
from Kleisli coinduction leads to observations and structures that are completely standard in the 
field. 

3.2 Markov Chains 

Our second example application, although conceptually rather a little more advanced, can be 
presented with much less technical detail, because we can build on standard accounts in terms of 
far more expressive mathematics; see for instance [9]. 

Consider discrete distributions on a set A, with possibly countably infinite support, TX = 
{ 71 - : A' —>■ [0,1] | n ( x ) = !}■ The distribution space TX has a convex structure. Hence its 

elements can be understood as formal sums J2iPi x i, and the operations 

Tf(Y,Pi x i) = I2Pif( x i) h( x ) = = E PiQij x ij 

i i i j ij 

define a monad T = (T, 77 ,/r) [10]. Alternatively, the elements n £ TX can be understood as 
stochastic row vectors n with columns index by A. By extension, Kleisli morphisms of type 
/ : A —> TY can be understood as right stochastic matrices F of shape A x Y, by currying. 

Lemma 3.11. Kleisli extension and composition are right matrix multiplication: 

r(n) = HF go Kl Wf = FG 

The T-coalgebras (A, e) are exactly the space-discrete, time-homogeneous Markov chains, a 
model class pervading almost all corners of science. The T-algebras are convex sets: structures (Y. f) 
with a set Y closed under a given interpretation / of formal convex combinations. In a very abstract 
but logically precise sense, the corresponding ca-homomorphisms employ convex sets as models 
of the long-term behaviour of a Markov chain: Let Y be a set of possible long-term behaviours 
(refraining from jumping to the conclusion of well-known candidates). A ca-homomorphism is a map 
h : X —> Y of Markov states to long-term behaviours that is consistent with transition structure e, 
namely e(x) = Y^iPi x i implies h(x) = f (Y^iPi^( x i))■ That is, the long-term behaviour associated 
with a state is an appropriately weighted convex combination of the long-term behaviours of its 
successors. 

The default candidates for long-term behavior that we have alluded to are of course the stationary 
distributions, the fixpoints of /*. That our ca-homomorphisms are a generalization, is apparent from 
the Kleisli-coinductive solutions: the matrix version by Lemma 3.11 of the coinduction property 
gives the linear fixpoint equation E' = E^E, such that the rows of are logically independent 
unitary left eigenvectors of E, which is the standard characterisation of stationary distributions. 


Note that, even though the intuition of distributions represented as formal convex combinations 
of states is quite dominant, the target space Y is completely abstract in this framework; any 
convex set will do from an axiomatic perspective. The classical presentation can be retrieved as the 
distinguished target T-algebra (TX,fix)- 

Regarding the degrees of freedom in ca-homomorphisms, Markov chains are an interesting 
escalation over exact dynamical systems: On the one hand, they support an analogous concept of 
coproduct structure, namely the so-called communicating classes, which are the strongly connected 
components of the non-zero transition probability graph. On the other hand, they have a well- 
known additional degree of freedom: even irreducible (totally connected) Markov chains can have 
ambiguous stationary distribution, if they are periodic; that is, some state has zero transition 
probability to itself in n steps, for arbitrarily large n. 

Note that, for subtle interpretative reasons, the direction of time in this subsection is reversed 
with respect to the previous one; the language of equation coalgebras speaks about the future rather 
than the past, whereas the algebra constructs the past rather than the future; in philosophical 
terms, the latter is governed by final rather than causal reasons. 


4 Applications of Co-Kleisli Induction 

4.1 Economic Games 

As the first application of co-Kleisli induction, we review the coalgebraic presentation of economic 
game theory from [1]. They construct perfect information games in extensional form by means of 
anamorphisms into the final coalgebra of a suitable game functor. We adapt their formalization 
slightly. 

Definition 4.1 (Game Functor). Fix a set A of agents and a set C of choices. The game functor 
is defined as follows: 

GX = R- 4 + A x X c 

Elements are of either of the forms b\(u : A —> R) or 02 (a,m : C —> X). They are understood as 
stating, respectively, that ‘the game terminates with payoff u(x) for all agents x £ A\ or that ‘the 
game continues with agent a’s turn and continuation m(c) for any possible choice c £ C that a can 
make’. 

In [1], possibly infinite game trees are unfolded from coalgebras (X, e) that model games as 
transition systems. The concrete states from X are erased by the anamorphism as usual, only the 
transition structure is retained. By contrast, we shall consider the cofree comonad G°°, and the 
lifted G ,00 -coalgebras (X, e^). In comparison, they can be understood as running ‘with logging’; the 
original states are retained as node labels in the tree. This allows for concise treatment of game 
tree evaluations as node label processors. 

We shall present only a single, very well-known evaluation algebra, in order to demonstrate the 
applicability of our approach. The field of game theory is broad and deep, and in great need of 
formal frameworks that give concise and elegant notations. Thus we mark the extension of the 
present sketch as an interesting topic for future research, which however requires a great deal more 
effort and space. 

One of the historically and logically most basic evaluation techniques for games is backward 
induction [ 22 ]. 

The typical application of backward induction, especially in economic uses of game theory such 
as discussion in [ 1 ], is the computation of expected payoff, under the assumption that all agents 
choose such as to maximize their own payoff. This notion can be stated extremely concisely as the 
generating G-algebra (R- 4 , k 0 ) of a cobasic G°°-algebra (R -4 , k = k 0 ° Xr-a): 

fco((-i(w)) = u fco(t 2 (a,m)) = m(argmaxc.m(c)(a)) 

It follows that any co-Kleisli-inductive solution Ad : G°°X —> R ' 4 satisfies Ad (< 7 ) = u if 0 x (g) = 
Li(u), and the fairly convoluted (k$ o m) (arg max c. (Ad o m)(c)(a)) if 9x{g ) = 12 ( 0 , rn), which is on 
second thought a natural formalization of the above prose description for game trees in general, 
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regardless of the game being played. The canonical extension to a ca-homomorphism from a 
concrete game coalgebra ( X , /), namely fc* o /, restricts the evaluation to legal trees with a given 
starting state. 

It is well-known and appreciated in game theory that there is a unique solution, if only finite 
games are considered. In our framework, that is to say that finite games can be accomodated by 
the recursive subcomonad of G°° [21]. We conjecture that some standard techniques for forcing 
canonical solutions in infinite cases, such as discounting , can be expressed conveniently in terms of 
a non-cobasic evaluation algebra. 

As a sidenote, the very term backward induction embodies the characteristic difficulty of writing 
about scientific modelling, namely the required double reading with regard to formal representation 
and phenomenological interpretation, respectively. As we have noted in the preceding section, the 
direction of time, and hence the word ‘backward’ is a concept of the phenomenological dimension; by 
contrast the algebraic structure and hence the word ‘induction’ is a concept of the formal dimension, 
making the composite term essentially a category mistake. 

4.2 Lindenmayer Fractals 

The following, last example application, while another instance of co-Kleisli-induction for a fixed 
algebra, has added value in two respects: Firstly, it demonstrates the definition of a domain-specific 
comonad as part of the formal modelling framework. Secondly and paradoxically, it makes perfectly 
reasonable use of essentially non-well-founded induction. 

It is a well-known statement, and frequent subject of mathematical instruction, that certain 
fractal shapes can be specified by Lindenmayer systems (L-systems), a formalism akin to Chomsky 
grammars that has originally been developed for formal modelling of the growth behavior of 
plants [11, 13, 15]. This statement, while morally true, is to be taken with a pinch of salt; 
considerable effort is required for a satisfactorily rigorous explication. We shall demonstrate that 
the framework of ca-homomorphisms is a natural and expressive background for this line of reasoning. 
Confer [7] for another account of fractals that builds on ca-homomorphisms (their diagram 3), but 
proceeds in an interestingly different way. 

Definition 4.2 (Curve). Fix some real Banach space V, say without loss of generality V = M 2 . 
Fix a real interval I. An I -curve is a continuous map /:/—>• V. 

Definition 4.3 (Unit Curve). Fix some unit vector e £ V with e ^ 0, say without loss of generality 
e = (1,0). An e-unit curve is a [0, l]-curve / with /(0) = 0 and /(1) = e. We write C\ for the set 
of all unit curves. 

Example 4.4. The simplest unit curve is the stroke: str(x) = xe. 

Definition 4.5 (Context-Free L-System). As observed in [20], simple context-free L-systems, with 
nonterminal symbols only, are finite coalgebras of the (nonempty) list functor CqX = X + : The 
carrier is the set of nonterminals, and the operation is simply the set of production rules. The 
traditional semantics of an L-system ( X , /) is the sequence of iterations (/*’ n (s)), seeded with a 
start symbol s £ X. 

We write () for the empty list, and x:: r for a list with first element x and rest list r, as in the 
programming language ML. Complex lists X\ x n :: () are condensed to (x \,..., x n ). The 

operator © denotes binary concatenation of lists. 

Definition 4.6 (Turn). Fix some group (G, •, i) of linear isometries on V, here the special orthogonal 
group of rotations p a = _ ™“) parameterized by angle a. The elements of G are called turns. 

The action on vectors is written • : G x V —> V. 

Definition 4.7 (Fractal L-System). L-systems for specifying fractals are extensions of simple 
context-free L-systems, that is, finite coalgebras of a more complex functor. In particular, they come 
with an additional real-valued shrink factor as a per-rule attribute, as well as terminal symbols for 
turns. As hinted in the short paper [20], these can be expressed conveniently in terms of composition 
£f r = SCqT with additional functors S = JR. x Id and T = G + Id, respectively. Note that £f r is 
monotonic and even preserves inclusions. 
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Example 4.8 (Koch Curve). The following singleton L-system specifies a fractal shape known as 
the Koch curve. We use a visually evocative, short notation where the shrink factor appears as 
superscript to the production rule arrow, and nonterminals and turns, represented by their angles, 
are simply juxtaposed. 

K -^4 K (+a) K (-2a) K (+a) K 

This contains the same information as the more formal, but also much more cumbersome full 
notation of / : {K} —> Cf r {K} as 

f(K) = (3, (i 2 (K), ^(r+a), i 2 (K), n(r- 2 a ), iz(K), n(r +a ), i 2 (K))) 

In order for an L-system to meaningfully specify a fractal curve, additional constraints must be 
met. These are not readily expressed in syntactic terms, and because of the inherently self-referential 
nature of L-systems, great care must be exercised not to state them in a circular way. 

Definition 4.9 (Well-Formed Rule). We give an abstract geometric interpretation of expressions 
of type £f r X, under the assumption that each nonterminal represents a unit curve. To this end, 
we track the overall change of position and orientation effected by sequentially tracing curves and 
executing turns, as in ‘turtle’ graphics, disregarding the shrink factor. 

span : CqTX V 

f 5 

span(t ) = < p • span(u) 

I e + span(u) 

A rule is called well-formed, if and only if it is properly shrinking, spans the unit vector when 
shrunk, and has balanced turns: wfr((a,t)) a > 1 A span(t) = ae A dir(t) = i. 

Lemma 4.10. The abstract interpretation is additive. 

span(t ® u) = span(t) + dir(t) • span(u) dir(t ® u) = dir(t) • dir(u) 

Lemma 4.11. Well-formedness is additive: wfr((a,t)') A wfr((b,u)) =4> wfr((a + b, t ® u )). 

Lemma 4.12. Well-formedness is invariant under the fractal Lindenmayer functor. 

This result allows us to abstract from the contents of a rule, retaining only the shape, and to 
construct a subfunctor for well-formed rules. 


if t = {) 
if t = L\{p) :: u 
if t = i 2 {x) :: u 


dir : C 0 TX — G 

U if t = {) 

dir(t ) = < p ■ dir(u) if t = L\{p) :: u 
I dir(u) if t= i 2 {x) :: u 


Definition 4.13 (Rule Shape). Let 1 = {*} be the final object, with unique morpliisms lx ■ X —> 1. 
We call £f r l the shape space, and £f r lx ■ £frX —$■ £f r l the shape map. 


Lemma 4.14. A ride is well-formed if and only if its shape is. 

Definition 4.15 (Well-Formed Fractal L-System). We obtain a subfunctor £ w f r = C{ T \ w f r of the 
general fractal Lindenmayer functor by restriction to we 11-formed rules. 


Definition 4.16 (Rule Interpretation). The concrete interpretation of rules differs from the abstract 
interpretation with span by actually composing the subcurves listed in a rule, concatenating their 
domains. We model this idea as a £f r -algebra (Y, draw) over the general function space Y = Y[°’ 1 ] ) 
not assuming continuity for now. 

In order to obtain the shape of a unit curve, both the domain and the result need to be 
compressed, according to the length of the rule and the shrink factor, respectively. In the following, 
let t denote the length of list t € CqTY. 


draw((a, t )) (2) 


stepff) {x) 


1 1 step(t)(£z) 

'0 

0 

fix) 

e 

p • step(u)(z — 1) 
e + step(u){z — 1) 


where 

if 2 € [0,1] t 

if 2 = 0 t 

if 2 G (0,1) t 

if 2 = 1 t 

if 2 > 1 t 

if 2 > 1 t 


A iP) 

u 

(A) 

G.if) 

: u 

(B) 

(/) 

: u 

(C) 

G if) 

: u 

(D) 

Gip) 

u 

(E) 

<-2 (/) 

: u 

(F) 
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The auxiliary function step has dependent type; for every argument t of length I > 0, it yields a 
function of type . Note that the case t = ( ) does not occur, even recursively. 

By overlapping the subcurves at their glueing points in the composition, and even sending turns 
to constant points, we depart radically from the way of [7], where injectivity is a major theoretical 
goal, and successfully solved technical obstacle. Whereas they have been concerned with faithful 
representation of fractal space, our approach is more a computational account: By virtue of the 
piecewise constancy introduced by turns, a point f(x) for random x can be computed exactly by 
finite recursion almost surely , except for some pathological L-systems. We discuss this property 
more deeply in a forthcoming companion paper. 

Lemma 4.17. The endpoint of step is predicted by span: step(t)(£) = span(t). 

Theorem 4.18. Interpretation of well-formed rules preserves unit curves: drawer) £ C\ for all 
r £ T w f r (C*i). 

Proof (Sketch). By piecewise continuity with agreement at glueing points. □ 

Definition 4.19 (Well-Formed Rule Interpretation). We obtain an £ w f r -algebra (Ci, draw w { r ) by 
restricting rule interpretation draw to unit curves. 

Definition 4.20. Let a : £( w )f r => V be the obvious natural transformation that extracts the set 
of nonterminal symbols occurring in a rule. Let P = p) be the well-known finite powerset 

monad. Now consider the cofree comonad over £ w f r . The object C^ fT X can be considered as the 
space of non-well-founded X-labelled trees of nested well-formed rules. From these we can extract 
sets of nestedly occurring rule shapes: The natural transformation 

• p = r] ox- ^wfr extracts the single root shape; 

• a = aF°° o 8 : £“ fr => extracts the immediate subtrees. 

Thus the extraction of all shapes occurring at some finite depth in t £ £“ fr X is defined as 
Sx(t) = DZoi'PJ x ° Px ° 0x")(*)) where Kleisli extension is over P. Confer the relation (~-> e ) 
from Definition 3.2. 

Definition 4.21. We obtain the shape-finitary subcomonad D of the cofree comonad of well-formed 
L-systems: DX = {t £ £“ fr X \ Sx{t) finite}. Clearly, all actual such L-systems (X,f), which 
are finite £ w f r -coalgebras, extend to D-coalgebras (X,/*). Note the analogy to rational coalgebras 
[2, 17]. 

Finally we construct an interpretation (coequation) algebra, the very goal of this whole subsection, 
by extending draw w f r to a cobasic £“ fr -algebra, and restricting its domain to shape-finitary trees: 
k = draw wil o \c 1 |dCi • 

The ca-homomorphisms h : (X,/) —> (Ci,fc) map nonterminal symbols of L-systems to unit 
curves, with the obvious self-similar geometric consistency condition. Now we can make the point 
of our nontrivial efforts and restrictions. 

Theorem 4.22. The coequation D-algebra ( C\,k ) is corecursive; a unique ca-homomorphism exists 
from any D -coalgebra. 

Proof (Sketch). By uniform continuity and uniqueness on the dense subdomain of recursive glueing 
points. □ 

Note that, while the comonadic structure appears not to be used by the cobasic algebra, it is 
required for expressing the restriction to shape-finitary L-systems, which in turn is a necessary 
condition for our proof; infinitary counterexamples can be conceived. 
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A Omitted Proofs 


This section gives or completes all proofs that we consider informative or non-elementary. Some 
boring immediate or inductive proofs are omitted altogether. 

Proof of Lemma 1.1. Let F preserve coproducts. Let index i range over some set / in all of the 
following. Then, for any family of objects (X,), there is a unique isomorphism j : FA', —> F JJ • X, 

such that, for any family of morphisms (hi : X, —► Y), we have [Fht] = F[hi\ o j. Now let (X*, ff) 
be a family of F-coalgebras. Then (X, /) = (]J- X i; j o ]J fi ) is again an F-coalgebra. Let (Y. g) 
be a fixed f-algebra. The ca-homomorphisms of type h : (A', /) — > ( Y , g) are in one-to-one 
correspondence with the families of morphisms (hi) where all members are ca-homomorphisms 
hi ■■ (Xi,fi) (Y,g), namely h = [hi]: 

II -—-> IL FX i ~^ F II Xi 



For products of algebras dually. □ 

Proof of Lemma 1.2 & 1.3. By natural laws of coproduct and monad, and catamorphic g. 

FF* 



For 0 dually. □ 

Proof of Lemma 1.4 & 1.6. The former clause by compatibility of / with /i; the latter clause by 
naturality of g. 


TTX TX 

Tf f 

TX -> X 

/ 


TTX 

TTh 

TTY 


rx 


r-Y 


y TX 

Th 

-t T Y 


For v dually. 

Proof of Lemma 2.2 & 2.11. Proposition laid out in underlying category. 


TXe 

Te t 

TTY 


(e f r 


X 


Ml 


-4T7 


O 


For kf dually. 


□ 
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Proof of Lemma 2.3 & 2.12. By composition of Lemma 1.4, second clause, with Lemma 2.2. 


■X 


TX 4— 

Te + (e+)^ e + 

tty ry 


TTh 


Th 


TTZ - y TZ 

For kf dually. 

Proof of Theorem 2.4 & 2.13. By composition of Lemma 1.4, first clause, with Lemma 2.2. 

TX^ -X 

Te t (eB* e+ 

TTY>TY 


□ 


Tf 

TY 


->y 


For Ad dually. 

Proof of Lemma 2.8 & 2.17. By naturality of k. 


n 



For x dually. □ 

Proof of Lemma 3-4- Assume that x e (t,x') and x (u,x r ). Then there are m,n such that 
e*’ m (x) = (■ t,x') and e* ,n (x ) = (u, x'). Assume without loss of generality m < n. Then 

(u, x') = e* ,n (x ) 

= ((e*) n ~ m o e*’ m )(x) 

= (e*r- m ((t,x')) 

= (t + v, x') where e*’ n ~ m (x') = (■ v , x') 

x' -+ e ( V,X') 


With v = 0 by consistency, we obtain t = u. O 

Proof of Theorem 3.7. By well-founded induction: Assume e'(x') has been chosen for all x' -< e x. 
Then either x is minimal, and e'(x) can be chosen arbitrarily; or otherwise there are some x',x",t 
such that x' -< e x, x" ~ e x, and e(x") = (t, x'). Then choose e*(x") = ( t+u , y) where e\x') = ( u , y). 
By cancellation and Lemma 3.5, necessarily e'{x") = e'(x) for all x" ~ e x. □ 
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Proof of Theorem 3.8. Consider the relation )~ e restricted to each equivalence class C of ~ e in 
turn, by Lemma 3.3. Either it is well-founded, such that Theorem 3.7 applies directly; or otherwise 
the choice of any element x partitions C into the down-set D of x, which is an infinitely descending 
chain, and its complement C\D, which is well-founded. Then fix u,y arbitrarily, and choose 
e'(x') = (u — t,y ) where by Lemma 3.4, t is uniquely determined by x ~-+ e (t,x r ), for each x! £ D. 
Proceed on C \ D as in Theorem 3.7. □ 

Proof of Theorem f.18. Consider a well-formed rule r = (a, t). 

The map step(t) is defined piecewise: glueing points at integer arguments are fixed, by cases (B) 
and (D). Consider the first pair of adjacent glueing points, at 0 and 1: they are either identical 
0, and the map is constant on the whole closed interval [0,1], by case (A), or they are 0 and e, 
respectively, and a unit curve is spliced continuously in between, by case (C). For further pairs, 
continuously shifted and turned analogs hold, by cases (E) and (F), respectively. Hence the map is 
continuous as a whole. 

Furthermore, we have step(t)( 0) = 0 by case (B), and step(t)(£) = ae by well-formedness and 
Lemma 4.17. Thus the rescaled map drawfr) is a unit curve. □ 

Proof of Theorem 4-22. By Lemma 2.8, the D-ca-homomorphisms h : (X,f$) —> ( C\,k ) can be 
understood as £ w f r -ca-homomorphisms h : ( X , /) —> (6\, drawer). That is, they are recursive in 
the nonterminal slots of a rule. Now proceed by considering the resulting curve pointwise, and case 
distinction in the workhorse function step: 

By the base cases (A), (B) and (D), the inductive definition of h in terms of induction step draw 
is well-founded, and hence the resulting curve defined uniquely, on a dense subdomain J C [0,1]; 
namely on isolated points of recursive interval subdivision by (B) and (D), and on closed subintervals 
by (A). By the induction hypothesis, the recursively included subcurves are unit curves and hence 
continuous on J. By finite amount of possible shapes, the extent of subcurves on J is bounded, 
but they are shrunk recursively by unbounded factors. Thus the resulting curve is totally bounded 
and hence uniformly continuous on J. By a standard result, a uniformly continuous map on a 
dense subdomain extends uniquely to a continuous map on its closure, the whole unit domain [0,1]. 
Clearly, that extension is a unit curve. 0 
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B Supplementary Figures 

This appendix illustrates the construction of fractal curves from L-Systems. Two example systems 
for very well-known fractals are considered. 

• The Koch curve 

K K IKrr K11\ 


• The Sierpinski triangle 


UIDrUrDl 
D rUlDlUr 


where U yields the standard upright triangle, and D its downward mirror image. 

The fractals are visualized by equidistantly spaced samples from the domain interval. The unit 
vector e spans the baseline from left to right. Note that analytically exact coordinates can be given 
for each sampled point by recursive symbolic interpretation. Sparse points are connected by dashed 
lines as a visual aid; these are not valid interpolations of the fractal. Note also that there may be 
fewer visible points than samples due to non-injectivity. 

The figures have been generated by programs written by an author in Haskell and R. The source 
code is available on request. 
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Figure 1: Koch curve h(K), sampled every 1/666 



Figure 2: Koch curve h(I \), sampled every 1/400000 
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Figure 3: Sierpinski triangle h(U), sampled every 1/666 



Figure 4: Sierpinski triangle h(U), sampled every 1/400000 
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